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ABSTRACT 

O 
O 

CN ' We present Chandra observations of EMSS 1358+6245, a relaxed cooling 

^ ■ flow cluster of galaxies at z = 0.328. We employ a new deprojection technique to 

pLn ' construct temperature, gas, and dark matter profiles. We confirm the presence of 

Q\ • cool gas in the cluster core, and our deprojected temperature profile for the hot 

component is isothermal over 30 kpc < r < 0.8 Mpc. Fitting the mass profile to 
CN ■ an NFW model yields = 153lg3^ ^P^ ^ ~ 8-4-2 s- We find good agreement 

^ ■ between our dark matter profile and weak gravitational lensing measurements. 

We place an upper limit of 42 kpc (90% confidence limit) on the size of any 
Q\ • constant density core. We compare this result to recent simulations and place a 

conservative upper limit on the dark matter particle scattering cross section of 
0.1 cm^ g~^. This limit implies that the cross-section must be velocity dependent 
• if the relatively shallow core mass profiles of dwarf galaxies are a direct result of 

dark matter self-interaction. 

O 

ti 

^ ■ Subject headings: X-rays : galaxies: clusters — cosmology : dark matter 

> 

X 

\-i ' 1. Introduction 

Galaxy clusters are powerful laboratories wherein to test the structure formation the- 
ories and simulations of modern cosmology. For example, while enjoying many successes 
in explaining a large number of observational results (Navarro, Frenk & White 1996, 1997; 
Moore et al. 1999), the cold dark matter (CDM) paradigm appears to be inconsistent with 
details of the structure of cluster dark matter halos (Spergel & Steinhardt 2000) and galactic 
rotation and density profiles (Spergel & Steinhardt 2000; Moore et al. 1999). The question of 
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whether a CDM universe can produce the observed mass profiles of galaxy clusters without 
serious modification is currently unresolved, although the emerging consensus is that CDM 
alone probably fails to produce cluster cores like those observed (e.g. Spergel & Steinhardt 
(2000); Dave et al. (2001); Yoshida et al. (2000); see Taylor & Navarro (2001), however). 

The launch of the Chandra X-ray Observatory (Weisskopf et al. 2000) has made it pos- 
sible to study the mass distribution in the central regions of galaxy clusters where the dis- 
crepancies between CDM models and observations are the most glaring (Moore et al. 1999; 
Firmani ct al. 2000). The ~ 0."5 spatial resolution of the Chandra/ ACIS-S instrument, 
which corresponds to about 3 kpc at z = 0.3, is now comparable to that of ground-based 
optical measurements, (although in practice low-N photon statistics limits the actual reso- 
lution achieved). For nearby clusters, for sufficiently long exposures, one can measure the 
mass profile of a cluster in the inner few kpc (David et al. 2001), scales which are generally 
inaccessible to gravitational lensing studies. 

The advantage possessed by gravitational lensing studies, however, is that the mea- 
surement is independent of the dynamical state of the gravitating matter. To derive a total 
gravitating mass profile from the baryonic X-ray emission map one assumes hydrostatic equi- 
librium (Sarazin 1988; Allen & Fabian 1994), which in practice means that one assumes that 
the cluster is supported via isotropic thermal pressure (i.e. frot/c <S 1 and |Bp/nA;T <^ 1), 
and that no recent merger event has caused a disruption in the pressure, temperature and 
density profiles. One would expect clusters with fairly circular isophotes, and which contain 
cooling flows (which are thought to be disrupted in merging events - see Allen & Fabian 
(1997), Allen (1998) and Fabian et al. (2001)), to merit the simplifying hydrostatic treat- 
ment. 

EMSS 1358+6245 (a.k.a. CL 1358+6245, ZwCl 6249) is such a galaxy cluster. Discov- 
ered optically by Zwicky (Zwicky & Herzog 1968) and later in X-rays with the Einstein 
IPC (Luppino et al. 1991), it is a nearly circular, apparently relaxed cluster, with an Abell 
richness class 4 (Luppino et al. 1991; Hoekstra et al. 1998) and a central bright cooling flow 
(Bautz 1997; Allen 1998). Because of its moderate redshift {z — 0.328), the entire cluster 
fits on the ACIS S3 chip, leaving enough blank field to derive a rehable background esti- 
mate. Figure 1 shows a broad band (0.3-7.0 keV) ACIS S3 image of EMSS 1358+6245, 
adaptively Gaussian smoothed on scales of 1 to 10 pixels (~0.5 to 5"). The contour levels 
are listed from inner- to outermost; contours associated with point sources in the field have 
been omitted. The presence of the cooling flow and the regularity of the X-ray isophotes 
suggest that hydrostatic equilibrium is a reasonable approximation of the dynamical state 
of the cluster gas. Figure 2 shows a mosaic of archived HST exposures of the central region 
of EMMS 1358+6245 with adaptively-smoothed soft-band (0.3-2.0 keV) X-ray isophotes su- 
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perimposed. Note that the peak of the emission (i.e. the central green dot), at a value of 
2.95 photons/(0."49)^, appears to be slightly offset from the position of the central dominant 
galaxy (CDG). We shall return to this point in §2. 

The mass of EMSS 1358+6245 has been measured using both weak (Hoekstra et al. 
1998) and strong (Pranx et al. 1997; Allen 1998) gravitational lensing, and estimated from 
ASCA and i? 0574 T observations (Baiitz 1997; Allen 1998). These measurements arc roughly 
consistent, with a mass of about 4 x lO"*^"^ Mq enclosed within the 1 Mpc. In this study, we 
use a Chandra ACIS-S imaging spectroscopy to derive a mass profile of EMSS 1358+6245 
and compare it to lensing measurements. 

We briefly outline the Chandra observation and data analysis procedure in §2. We de- 
scribe our deprojection and modelling technique in §3, and compare our results to previously 
published measurements in §4. We then use our profile to constrain possible dark matter 
matter candidates in §5. Finally we summarize our findings in §6. We assume flm = 1, 
Qa = 0, and Hq = 50/i5o km 8"^ Mpc"^ throughout (1" = 5.79 kpc at z = 0.328 for h^o = !)• 

2. Chandra Observations, Astrometry, and Data Analysis 

EMSS 1358+62 was observed with the Chandra X-ray Observatory for approximately 
55 ks on 3-4 September 2000 using the S3 chip on the ACIS detector. The level 2-processed 
data were aspect-corrected and filtered for periods of high background using the CIAO soft- 
ware package, according to the standard procedures described in the CIAO Science Threads 
(2001). 

The center of the projected X-ray emissivity (Figure 2) is offset from the optical position 
of the CDG by about 2" (about 12 kpc). The offset, if real, has important consequences 
for the dynamical state of the cluster, so we carefully examined the the astrometry of the 
Chandra field. We ran the CIAO wavelet source detection routine wavdetect (CIAO Science 
Threads 2001) on the ACIS S2/S3 fields (i.e. those chips closest to the ACIS-S aimpoint), 
and cross-correlated the output with the Tycho-2 (H0g et al. 2000) and USNO-A2.0 (Monet 
et al. 1998) catalogues. Five USNO sources and zero Tycho sources were within 3" of a 
wavlet-detected Chandra source. Figure 3 shows the location of the five USNO sources in 
relation to the Chandra ACIS S2/S3 field, the HST mosaic, and the CDG. Although the 
Tycho-2 catalogue is nominally more accurate than the USNO-A2.0 catalogue, with internal 
astrometric standard errors in general less than 90 mas (H0g et al. 2000), above S = —20 
the USNO positions are usually within 20 mas of the Tycho-2 positions (Assafin et al. 2001), 
making USNO-A2.0 as effective an astrometric standard. Of these five sources, two appear 
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somewhat extended, either because they are galaxies or because they are sufficiently far from 
the ACTS aim point. The remaining three point sources were used as astrometric anchors 
to calculate the difference between the Chandra and USNO coordinates. These offsets are 
shown as filled circles in Figure 4. The mean offset is represented by the small open circle, and 
its error by the ellipse. About half of the CDG/X-ray peak offset is due to the Chandra field 
astrometry. Figure 5 shows the center of the cluster including this astrometric correction. 
The top panel shows broad-band (0.3-7.0 keV) X-ray contours superimposed on the 0-10 
pixel adaptively smoothed, 2 x 2-binned broad-band Chandra image. The same contours 
are also shown on the unsmoothed 2 x 2-binned Chandra image (middle) and on the central 
field of the HST mosaic (bottom). The offset between the CDG and the X-ray peak that 
remains after the astrometric correction is 0."99, while the error in the correction is 0."97. 
Therefore, the CDG and the X-ray emission peak are cospatial at about the 1-cr level. 

After locating the center of the projected emissivity we divided the cluster into 10 
concentric annuli such that each contained a minimum of 2000 counts (1600-2000 after 
background subtraction), the center annulus being a disk. Owing to decreased signal-to- 
noise in the outer bins, annuli 9 and 10 were constructed to contain roughly 4000 (2700) and 
10,000 (3300) photons, respectively. We found that the extra counts in the last two annuli 
were required for numerical stability during the iterative fitting procedure, and they also 
helped to reduce the "noise" in the derived temperature profile. The annuli ranged (outer 
radius) from 11.5 to 285.8 pixels, which translates to 5.66 to 140.56 arcsec. At z — 0.328, 
this corresponds to ~ 30 to 800 kpc, for Hq — 50 km s"^ Mpc~^ and Qq — 0.5 (see Table 1). 

All obvious point sources were removed, and a spectrum was extracted from each an- 
nulus. RMF and ARF response matrices were constructed according to the CIAO Science 
Threads (2001), as was a background spectrum from an annulus exterior to the outermost 
bin. The spectra were recorded in PI format, and grouped such that there were a minimum 
of 20 counts per channel. 



We assume that the cluster is a spherical, self-gravitating pressure-supported plasma 
whose X-ray emission is optically thin. We assume hydrostatic equilibrium throughout (ex- 
cept for small blobs of cold gas in the innermost regions which, although not hydrostatic, are 
assumed to be in a steady state). The hydrostatic equation can be written (Sarazin 1988) 



3. Cluster Deprojection and Mass Determination 
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where T and p are the local (baryonic) gas temperature and density, r is the spherical radius, 
M(r) is the total mass enclosed within r (i.e. baryons plus dark matter), and rUp and are 
the proton mass and mean particle weight, respectively. 

To check our assumption of spherical symmetry we examined the surface brightness 
map for deviations from circular symmetry. Specifically, we performed a simple analysis to 
estimate the ellipticity e (and position angle 9) of the cluster elongation. Wc calculate the 
luminosity of a pair of wedges A, centered on the X-ray peak of the cluster, and compare it 
to a second pair B perpendicular to the first (see Figure 6). The wedges have semiwidth 5, 
radius Ra (or Rb)-, and phase (or 6' + 7r/2). The luminosity of the two wedge pairs is then 
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Maximizing La{9) for Ra/ Rb = 1 yields the position angle, and then scaling Rb/Ra until 
La/ Lb = 1 gives an estimate of the axial ratio, which is simply 1 + e. The noise in the 
surface brightness map that tends to inflate La/ Lb for some 9 is amehorated somewhat by 
the smoothing introduced by the finite window size 5. Figure 7 shows the axial ratio and 
position angle as a function of Ra- The fiattening is consistent with an ellipticity of 0.34 
over an order of magnitude in Ra, for 5 — 7r/36. Because much of this fiattening is due to 
substructure in the soft emission due to the cool component (which shows an axial ratio of 
2.5 for Ra = 200), we feel confident that the assumption of spherical symmetry, especially 
in the hot, pressure-dominating component, is not unreasonable. We shall, however, return 
to the implications of e = 0.3 in §5. 

Se also fit a 2-D f3 model to the cluster emission map using the Sherpa data analysis 
package (CI AO Science Threads 2001). The profile is described by 
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r = ^u^ + v'^/{l-eY , (5) 
u and V are the rotated coordinates: 



f u\ ( cos 6* + sm.6 \ f X 
\ V J \ — sin ^ + cos J \ y 

and e and 9 are the isophotal eUipticity and position angle, respectively. (We have assigned 
(x, y) = to the X-ray emission peak.) Here ro is the core radius and /o is the central 
surface brightness. We fit for the core radius, the amplitude, and the isophotal eUipticity 
and orientation. We use the maximum likelihood method of Cash (1979), as implemented 
in Sherpa (ClAO Science Threads 2001), rather than a minimization, because the vast 
majority of the image pixels contain fewer than ~ 20 photons. 

We find an eUipticity of 0.174±0.008 for the cluster (the position angle has an oblate/prolate 
degeneracy), although the uncertainty is underestimated due to parameter correlations. This 
indicates again that spherical symmetry is a reasonable approximation. The best fit param- 
eter values are tq = 130 ± 30 kpc and f3 = 0.69 ± 0.02. 

In order to derive s^/ienca/ radial profiles we construct a model consisting of concentric 
spherical shells whose inner and outer radii correspond to the inner and outer cylindrical radii 
of the projected annuli in the data set. (Hereafter annuli/shclls will be labeled 1 through 
N, in order of increasing radius.) Quantities can be mapped between the two geometries 
through the upper diagonal matrix V, whose elements Vij contain the volume of spherical 
shell j intersected by a cylindrical shell formed by the projection of annulus i. If the volume 
emissivity of the gas in spherical shell j is Uj, then the luminosity of annulus i, Hi, is simply 




H = V.u, (7) 

where we have written the relation for all annuli simultaneously. Since H is related to the 
surface brightness / through Hj = 27r(r|— r|_^)/, the deprojcction of the surface brightness is 
accomplished by inverting V, which is completely specified by the binning geometry. Because 
V is upper diagonal, the solution set of volume emissivities is easily obtained working from 
the outer- to the innermost annulus. The emissivity of the outermost shell is simply un = 
Hn/Vnn: while the rest of the shells are calculated moving inward using 
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N 




i=j+i 



(8) 



The volume emissivity of each spherical shell is characterized by a temperature and 
a normalization. The X-ray emission from each shell is modelled spectroscopically using 
MEKAL (Mewe, Gronenschild & van den Oord 1986; Mewe, Lemen & van den Oord 1986; 
Kaastra 1992; Liehdal, Osterheld & Goldstein 1995) model in the XSPEC software package 
(Arnaud 1996), which describes the emission from an optically thin thermal plasma. 

Since a strong coohng flow signature was seen in both ASCA and ROSAT observations 
(Bautz 1997; Allen 1998), we first sought to verify its existence in the Chandra data. We 
modelled the inner two annuli separately using the cooling flow model CFLOW of Mushotzky 
&; Szymkowiak (1988) in XSPEC, together with a MEKAL component and an intervening 
absorber. The cooling flow model is described by a mass deposition rate M, a hot gas 
temperature and a low-temperature cut-off. We assume that the flow cools from the ambient 
cluster gas, so we pinned the hot gas temperature to the MEKAL temperature during the 
fltting procedure. We also assume that the emissivity is proportional to the inverse cooling 
time at the local temperature. An acceptable fit (x^/d.o.f. = 62.37/69) was obtained for 
a mass deposition rate of 40.1^18^9 Mq y~^ (90% confidence interval) and a low-T cut-off 
~ 1 keV. This result differs significantly from that of AUen (1998), who finds M = 690l29o 
Mq/j. The discrepancy is due to the extra (internal) absorption component which that 
study assigns to the cooling flow emission. While we flnd no particularly compelling reason 
to incorporate an additional absorber in this type of model (Arabadjis & Bregman 2000), we 
do not argue this point in the present study, since we merely wish to establish the presence 
of a cool emission component in the Chandra data. 

The bulk of the emission in the flow originates from the coolest gas, and so in the 
interest of simplicity we modelled the cooling flow component in the complete cluster analysis 
as a second MEKAL component, rather than with the CFLOW model. To determine the 
appropriate spatial extent of the cool gas we tried adding a second emission component 
(at lower temperature) to several of the inner spherical shells. (We shall refer to these two 
temperatures as and Tc, for the hot cluster gas and the cooling flow gas, respectively.) 
We ran models with two components in (A) shell 1, (B) shells 1 and 2, (C) shells 1 through 
3, and (D) 1 through 4. Model B had a lower reduced chi-squared value than did model 
A, and the temperature and normalizations of cool components in models C and D were 
essentially unconstrained. We therefore adopted B as our working model, and note that 
this suggests that the cool gas extends out to about 70 kpc. It should also be noted that. 
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unlike other deprojection methods, this method makes no assumptions about the shape of 
the temperature and density profiles other than spherical symmetry and the finite width of 
radial bins. 

The number of concentric shells that can used to construct this model is limited by two 
quantities, the memory allocation of the XSPEC program and the number of photons above 
the background in the spectrum. An active XSPEC model is limited to 1000 parameters, 
with only 100 allowed to vary at any time. Each MEKAL component contains 6 parameters, 
but we allow only kT and the normalization to vary. In addition, we include an intervening 
column of Galactic material which absorbs all emission components equally. For N concentric 
annuli in the data, there are N + N2 emission components, where N2 is the number of two- 
component shells. Including the absorbing column yields a model with 6(A'" + A^2) + 1 
parameters. This model is applied to N data sets, with the normalizations across them 
scaled by the appropriate geometric factors (i.e. ratios of V terms). Thus the complete 
XSPEC model contains A^[6(A^ + A^2) + 1] components, with 2(A^ + A^2) + 1 of them variable. 
In our case we had N2 — 2 two-component shells, thus limiting N to 11. 

In practice, however, we found that we had to limit N to 10 because the finite number 
of photons in our 55 ks spectrum caused instability during the numerical iteration if the 
temperature of any given shell was sufficiently unconstrained. It turns out that > 1800 
source photons per annulus are required to bring the "noise" in the fitted temperature profile 
to a level where numerical stability is achieved throughout the iterative process, which in 
our case resulted in 10 annuli (using A^2 = 2). Thus our model contains a total of 730 
parameters, with 25 of these floating, well under the XSPEC limit. 

The best-flt model has a chi-squared of 942.4 for 891 degrees of freedom (there are 916 
spectral bins in the data). The temperature and density profiles arc shown in Figure 8. The 
model shown at left allows the temperature of each shell to vary. Note that the temperature 
profile is consistent with isothermality, with an error- weighted average temperature of 7.16± 
0.10 keV, while the right side pins all the shells at 7 keV. In each case the second emission 
component is shown as the pair of points near R— 0.1 and 0.2 arcsec. The inferred density 
profiles are nearly indistinguishable. 

The temperature and density profiles can be used to derive an entropy profile of the 
cluster. The discovery of an "entropy floor" in galaxy clusters (Ponman, Cannon & Navarro 

1999) and its subsequent measurement in many clusters (Lloyd-Davies, Ponman, & Cannon 

2000) have led groups to attempt to explain it in terms of an additional source of heat 
(SNe and/or AGN activity) which acts either to preheat the gas that collapses to form the 
cluster, or which heats the cluster throughout its formation (Balogh, Babul & Patton 1999; 
Loewenstein 2000; Wu, Fabian & Nulsen 2000). The speciflc entropy of the cluster baryons 
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can be calculated using 



S = Tnf/^ (9) 



where Ue is the electron number density of the gas. In Figure 9 we plot the entropy profile 
of the cluster. The cold gas entropy is shown for the inner two shells as well. The entropy 
of the hot cluster gas toward the center is consistent with the entropy floor of 70-140 keV 
cm^ reported by Lloyd-Davies, Ponman, & Cannon (2000) in their sample of 20 clusters. 

The temperature and density profiles can also be used according to equation 1 to measure 
the total gravitating mass in the cluster. Figure 10 shows the mass profiles corresponding 
to the temperature/density profiles of Figure 8, with the logarithmic derivatives calculated 
as simple differences. The data points in the left panel represent M(r), the mass contained 
within spherical radius r, calculated using the density and unrestricted temperature profile 
shown in the left panels of Figure 8. The solid blue curve represents the best-fit "universal 
density profile" of Navarro, Frenk & White (1996, 1997), hereafter the NFW profile, a two- 
parameter family of models whose functional form is 

where is a scale length, po = ScPcru{z), and 6c is a characteristic (dimensionless) density 
dependent upon the formation epoch of the dark matter halo. Pcrit{z) is the critical density 
at the observed redshift; for a matter-dominated, Q^q — 1, h. — Q universe this is 

Pc..* = ^ (1 + ^)=^ (11) 

The total mass scales as so our baryon fraction is weakly dependent upon h. Of course, 
radial distances such as scale as . 

We report our best-fit values as and c, the concentration parameter, defined by 

^ ^ ^ ^ 200 

" Peru 3 ln(l + c) -c/(l + c) ^ ' 



In practice, because the form of our mass profile is M(r) and not p(r), we use the 
integrated mass of the NFW profile as our fitting template, rather than the density: 



M(r) = Mo [ln(l + r/r,) + (1 + r/r,)-^ - 1 ] 



(13) 



where Mq = ATipQr^. Using M(r) rather than p{r) spares us one numerical differentiation, 
which would augment the profile noise considerably. We find = ISS^gg^ Icpc and c = 

Kelson ct al. (2000) have measured the central velocity dispersion of the galaxy located 
at the cluster center using the Keck I 10m telescope. Using the low-resolution imaging spec- 
trograph with a square 1.05"aperture, they infer a central velocity dispersion of 311±11 km 
s~^. Prom HST photometry they obtain an effective radius of r = 3.9" (23 kpc). Asssuming 
a singular isothermal density profile with an istropic velocity dispersion, the enclosed mass 
is given by: 

M,so{r) = ^ (14) 

The dashed black line at left (Figure 10) represents an isothermal sphere at u = 311 km s~^, 
with the Kelson et al. (2000) measurement shown as a filled circle plotted at the effective 
radius. Although the point barely falls within the l-cr envelope, it should be remembered 
that this point is a model extrapolation from the aperture radius. If one uses the aperture 
radius rather that the effective radius, the point lies very near the NFW best-fit line. We 
shall return to this point in §5. 

The sohd blue curve at the right half of Figure 10 shows the best-fit {x^/d.o.f. ~ 2.1/6) 
spherical NFW profile projected along the line of sight, i.e. M{r) — > M[R) (in order to 
compare with the surface density maps derived from weak lensing measurements; see §4). 
Integrating equation 10 along the line of sight yields (Bartelmann 1996) 



{ln(a;/2) + :^7=f arctanhVl — x'^ x < 1 

ln(a;/2) + l x^l (15) 

ln(x/2) + yjf3j arctanVa;2 _ i x > 1 

where x — R/ts- 

Because of the "noise" in the solution, one point in the mass profile is unphysical. This 
is caused by a l-cr jog in the temperature profile at the seventh annulus. The (positive) 
logarithmic temperature gradient between annuli 6 and 7 swamps the density gradient (see 
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equation 1) causing M(r) < 0; this is represented in Figure 10 by an open dashed box at an 
arbitrary ordinate. 

In order to ascertain whether this jog represents an unseen feature in the data or a 
fundamental flaw in the method, rather than an ordinary statistical fluctuation, we varied 
the initial conditions in the fitting procedure to see if the iterative process would converge 
to a solution lacking such a temperature anomaly. After many attempts we found that we 
could only move the jog to a different annulus. As a more rigorous check on the method, we 
analyzed a cluster simulated with MARX that contained 10^ photons and no background (J. 
Houck & M. Wise, private communication). We derived a mass profile for this cluster using 
12 annuli, and found two unphysical points due to jogs in the otherwise smooth temperature 
profile. We therefore conclude that the single unphysical point in our mass profile for EMSS 
1358+62 is not due to undetected peculiarities in the data, but is more likely a nonexcep- 
tional statistical fluctuation. We note here that in their deprojection analysis of Hydra A, 
David et al. (2001), who used a smooth parametric temperature profile, were still forced to 
effectively smooth across radial bins in order to reduce the "noise" in the mass profile, even 
though they had roughly 5 times the number of photons at their disposal. 

Nonetheless, we ran the model again, pinning Th at 7 keV, to see what effect this would 
have on the mass profile. The resulting temperature and density profiles are shown on the 
right side of Figure 8. The effect on the mass profile, and its NFW fit, is small; the scale 
length increases from ~ 153 to 197 kpc, and concentration decreases from 8.4 to 7.0. The 
significance of these values will be discussed in §5 below. 

The gas mass and total gravitating mass profiles are shown in Figure 11, for the model 
wherein Th is allowed to vary. The Th = 7 keV model is very similar. The gas mass fraction 
rises from ~0.025 at 50 kpc to 0.18 at 0.6 Mpc, scaling as h^^ '^ (White & Fabian 1995). 
This is consistent with the study of White & Fabian (1995), who find baryon fractions of 
0.10-0.22 for a sample of 13 clusters (90% confidence). The mass of an isothermal sphere 
based on the Kelson et al. (2000) measurement, and its corresponding profile, is also shown, 
as a filled circle and a dashed line, respectively. Although these data do not allow us to 
distinguish between the gas, stars, and dark matter associated with the CDG, they clearly 
demonstrate that the total galaxy mass dominates the gas mass within the cluster center. 

Our model does not allow us to determine the morphology of the cool gas within the 
inner two shells of cluster gas, although from our ellipticity estimation we suspect that it 
is fairly aspherical. However, its mass, density, and volume are completely specified, since 
we have assumed pressure equilibrium between the two temperature components. For the 
unconstrained temperature profile model, roughly 10% of the core volume is filled with cold 
gas, dropping to a tenth this value at the adjacent shell. Table 2 lists the relevant quantities 
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for both gas components in the inner two shells. We note that the broadband (0.3-7.0 
keV) X-ray luminosity of the cool gas component in the inner shells is 4.4 x lO'^"' erg s~^, 
comparable to the luminosity of the "X-ray overluminous" ellipticals studied by Vikhlinin 
et al. (1999). 

4. Comparison with Lensing Studies 

In order to assess the robustness of our derived mass profiles, we compare them to the 
the weak lensing profile of Hoekstra et al. (1998). While gravitational lensing measurements 
may suffer from biases introduced by distant, uncorrelated large-scale structure, the effect 
on the mass determination is expected to be small for distant, rich clusters; for example, 
Hoekstra (2001) finds a typical la uncertainty of about 6% for measurements out to 1.5 h^Q 
Mpc. In addition, Allen (1998) has shown that lensing and X-ray mass determinations are 
consistent provided one takes into account the presence of the lower temperature gas in the 
coohng flow. 

Because weak lensing mass profiles measure mass within a projected radius, we integrate 
the best-fit NFW profile (equation 10) along the line of sight (equation 15). For both our 
treatments of the temperature profile of the baryons (i.e. variable vs. frozen at 7 keV), 
the resulting mass profile agrees well with the weak lensing profile. The strong lensing 
measurement (Pranx et al. 1997; Allen 1998) at 120 kpc exceeds the X-ray and weak lensing 
measurements by a factor of ~ 1.6. We note that although the lensing results we quote and 
our X-ray measurement assume spherical symmetry, Hoekstra et al. (1998) find evidence for 
an elliptical mass distribution with an axis ratio of ~ 0.3. Our analysis shows that while the 
X-ray surface brightness is very nearly circular at large radii, the axis ratio in the core could 
be as small as 0.7. The assumption of spherical symmetry leads to an average overstimate 
of a factor of 1.6 in strong- lensing analyses of simulated clusters, which are generally not 
spherically symmetric (Bartelmann 1995). 

Following Allen, Ettori & Fabian (2001) we compute an "effective velocity dispersion" 
of our best-fit mass profile. Using the NFW halo virial mass in the singular isothermal 
profile we obtain a = y/50H{z) VgC — 700l;g4^ km s~-^. This is substantially lower than the 
observed line-of-sight velocity dispersion of Fisher et al. (1998), 1027 ±50 km s~^. As noted 
by Fisher et al., the optical velocity dispersion is probably inflated by substructure along 
the line of sight to the cluster. As expected, our effective velocity dispersion agrees, within 
measurement errors, with the value of cr = 780 ±50 km s~^ derived by Hoekstra et al. (1998) 
by fitting an isothermal sphere to the weak lensing data. 
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5. Constraints on Self-interacting Dark-Matter 

Much has been made of the failure of CDM simulations to reproduce the observed 
structure of dark matter halos. A particular example of the problem is that CDM halos are 
steeper than observed dwarf galaxy halos. Parameterizing the core profile as p(r) = Po{r/rQ)" 
as r — > 0, CDM simulations suggest either a — —1 (Navarro, Prenk & White 1996, 1997) or 
—1.5 (Moore et al. 1999; Fukushige & Makino 2001). Ha observations low-surface brightness 
galaxies suggest a flatter profile, a = —0.5 (Swaters, Madore & Trewhella 2001; Dalcanton 
& Bernstein 2000). Instead of a triaxial, central cusp, as predicted by CDM simulations, 
Tyson et al. (1998) find a spherical = 35h^^ kpc core in the strong gravitational lens 
CL 0024+1654 using a multiply- imaged background galaxy (see Broadhurst et al. (2000) 
and Shapiro & Iliev (2000), however). A weak lensing study by Small et al. (2000) also 
indicates soft cores CL 1455-1-22 and CL 0016-1-16. The literature survey of Firmani et al. 

(2000) suggests that the central density of clusters is about 0.02 Mq pc"^, regardless of halo 
mass, whereas CDM predicts values of > 1 Mq pc~^ for dwarf galaxies, and larger still for 
more massive halos. Spergel & Steinhardt (2000), Firmani et al. (2000), Dave et al. (2001) 
and others use these discrepancies to argue for the existence of self-interacting dark matter 
(SIDM), which will tend to reduce both the core profile slope and central density of dark 
matter halos. 

X-ray measurements of cluster dark matter profiles, however, fall somewhat more in 
line with the CDM simulations. Tamura et al. (2001) find a central slope of ~ —1.5 for 
Abell 1060. Markevitch et al. (1999) find that NFW provides a good description of the mass 
profiles they derive for A2199 and A496 from ASCA and ROSAT observations, over the 
entire range for which they are able to derive a temperature profile. Allen, Ettori & Fabian 

(2001) observed A2390 using the ACIS S3 detector on Chandra and derive a mass profile 
which is consistent with the NFW profile, and find a concentration c — A, as predicted by 
CDM. Schmidt, Allen & Fabian (2001) also find adequate agreement between their mass 
profile of A1835, derived from a Chandra ACIS-S3 observation, and the NFW profile, and 
consistency with weak lensing measurements. They calculate a scale radius of 640^^20 ^P<^ 
and a concentration of 4.0lQ'g4. They find a slightly better fit when they model A1835 as a 
non-singular isothermal sphere, deriving a core size for the cluster (distinct from the scale 
radius of an NFW fit) of Vc = ^5^% kpc. David et al. (2001) determine a central slope of 
a ~ —1.3 for Hydra A, and find no evidence for a flat, SIDM dominated core, down to scales 
of ~ 40 h^Q kpc. However, their NFW fit to their mass profile results in a concentration 
c — 12, a factor of 3 larger than that predicted by CDM. 

In general we find that the NFW model provides an adequate fit to our mass profile. 
Like the David et al. (2001) study, however, our fit to the mass profile of EMSS 1358-1-6245 
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results in a concentration which is somewhat higher than that predicted by NFW for rich 
clusters. The appendix of Navarro, Frenk & White (1997) provides a prescription for calcu- 
lating cluster parameters at the collapse (i.e. formation) redshift Zcoii, which is a function of 
cluster mass in this model. Using the Einstein-de Sitter simulation with a CDM spectrum 
normalization ag = 0.63 and h = 0.5 (the SCDM model of Navarro, Frenk & White (1997)) 
we calculate a collapse redshift of Zcoii = 4.9; it is doubtful that a cluster this massive could 
have formed so early. NFW predicts a concentration of about 4 for massive clusters, but the 
scatter in their suite of simulations could be as high as 2. 

One might ask at what level SIDM can be ruled out by current observations. Recently, 
Yoshida et al. (2000) simulated cluster-sized halos and found that relatively small dark mat- 
ter cross-sections {(Tdm — 0.1 cm^ g~^) produce relatively large (40 kpc) cluster cores. 
Dave et al. (2001) simulated galaxy-sized halos and required interaction cross-sections as 
large as 5 cm^ g~^ to reproduce the shallow central slopes of galaxy density profiles. As 
noted by Dave et al. (2001) and explored by Firmani et al. (2001), these findings can be 
reconciled if the dark matter interaction cross-section is velocity-dependent. For example, 
the two findings are roughly consistent with 



with ao — 1 cm^ g~^, Vq — 100 km s~^, and a — —1. Hennawi & Ostriker (2001), however, 
use supermassive black hole demographics to rule out SIDM cross sections as large as this. 
They present a model with a = and do = 0.02 cm^ g~^ that is consistent with observed 
supermassive black hole masses. While their model does not remedy the galactic-scale dark 
matter halo problem, they sketch a scenario wherein the density cusps in such halos can be 
softened through black hole mergers. 

Our mass profile for EMSS 1358-1-6245 shows no evidence for a flat core on scales larger 
than 40 h^Q kpc. In Figure 12 we show power law fits to the mass profile as a function of 
the outermost point used in the fit. The profile is consistent with Mr oc r^'^, or a = —1.9, 
somewhat steeper than CDM predictions. In fact, the inner profile displays behavior which 
is the opposite of what would be required for nearing a fiat core; a should approach, not 
diverge from, as r — > 0. (The restriction is even more stringent since our model uses 
Hq — 50 km s~^ Mpc~^.) We place an upper hmit on the core size by modehng the profile 
as a softened isothermal sphere: 
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We find a core size Vc — 121^2 ^50^ kpc (90% confidence; /d.o.f. — 4.2/6). The contribution 
of the CDG to the mass profile for r > 30 kpc is sufficiently small (Figure 10) that we 
can directly compare our core size limit, r < 42/^50^ kpc, with Yoshida et al. (2000). The 
most conservative limit we can place must include the flattening that we determined in §3. 
Although the cluster as a whole exhibits an ellipticity of 0.34, it could be higher toward the 
core (see Figure 7. The worst-case scenario would be that our innermost annulus just grazes 
the minor axis of a core; in that case our limit should be scaled upward by a factor of 1 + e. 
Furthermore, in keeping with a conservative posture, if we assume that the axial ratio in the 
core (near r = 35 kpc; see Figure 7) is as high as 7/4, then our core fimit is 74 kpc, allowing 
us to conservatively rule out SIDM with Udm > 0.1 cm^ g~^. 

In order to strengthen the case for the core size upper limit, we verified that our depro- 
jection and modelling technique can indeed detect a core if one is present. We created a set 
of five simulated galaxy cluster observations to which we applied our deprojection and fitting 
method. Each simulated cluster follows a j3 model surface brightness profile, with (3 = 0.69. 
The (3 model core radii vq ranged from 6 to 96 ACIS S-3 pixels (34.7 to 556 kpc), increasing 
by a factor of 2 for each successive cluster. Otherwise the simulated clusters are similar to 
EMSS 1358+6245 - z = 0.328, and the photons are sampled from a distribution at T = 7 
keV. In addition, each of the 10 annuli contains ~ 2000 photons. Once we have fit for the 
temperature and gas density of each shell, the mass profile computed using the hydrostatic 
equation, which is then fit using a non-singular isothermal sphere. The resulting fits, shown 
in Figure 13, demonstrate that we can detect core structures in the mass profile. (The agree- 
ment between tq and Tc might at first seem too good, particularly when one considers that Tq 
is a projected beta model core radius, while Vc is a spherical softened isothermal sphere core 
radius. However, we are pinning down the shape of the mass profile near the core radius, 
where the non-singular isothermal sphere is actually a good approximation to a beta model. 
At r = Tc, for /3 = 0.69, the agreement is better than 5%.) 

In Figure 14 we plot the SIDM cross section constraints of several groups. The hatch 
marks delineate the region excluded by each constraint. The parameter space external to 
the triangular region is ruled out by Hennawi & Ostriker (2001) by requiring that SIDM 
(1) removes the halo cusp in dwarf galaxies, (2) does not cause core collapse in 20 km s~^ 
halos, and (3) does not cause galactic-mass halos to evaporate in a Hubble time. Hennawi 
& Ostriker (2001) rule out other areas not shown here by considering the formation of 
supermassive black holes (SMBHs), specifically the lack of a bulge/SMBH in M33, and the 
presence of a bulge/SMBH in the Milky Way. They arrive at SIDM which is unable to 
remedy the cuspy halo problem of dwarf galaxies, and and we direct the interested reader to 
that study. We have indicated the results of computer simulations of Yoshida et al. (2000) 
and Dave et al. (2001) (see arrows). Our mass profile of EMSS 1358-1-6245, in conjunction 
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with the Yoshida et al. (2000) simulations, provide the constraint (7d^(1000) < 0.1 cm^ g~^, 
ruhng out the space indicated in bhic. We observe that the core collapse and halo evaporation 
constraints of Hennawi & Ostriker (2001) appear to rule out the model of Dave et al. (2001) 
in which SIDM is responsible for dwarf galaxy halo cores. We note in passing that our 
observations are consistent with, though less restrictive than, the constraints derived by 
Hennawi & Ostriker (2001) from SMBH demographics. 

While our central profile is consistent with NFW, it is intriguing that the average density 
within the inner 50 kpc of our mass profile is ~ 0.025 Mq pc^^, consistent with the Firmani 
et al. (2000) sample. If the Kelson ct al. (2000) measurement of the (assumed isotropic) 
velocity dispersion at 0.55" is extrapolated to = 3.9" (22 kpc) as described in section §3, 
the resulting density is, perhaps coincidentally, p ~ 0.022 ± 0.005 Mq pc~^, suggesting that 
the data allow the density profile to be flat from 20 to 50 kpc. It is important to remember, 
however, that this extrapolation assumes that it is appropriate to characterize the CDG 
density profile as a singular isothermal sphere over the range 1.05" < R < 3.91". (More 
conservatively, if the isothermal model is assumed to hold only at the aperture radius, the 
profile is consistent with the absence of a core.) 

However uncertain, this result is interesting because it is contrary to what would be 
expected from a dark halo associated with the central galaxy, and may offer a clue to the 
existence of SIDM which flattens cluster cores on scales of 40 kpc. The issue can be addressed 
with a deeper Chandra observation to enable one to probe the mass profile in X-rays down to 
the effective radius of the Kelson et al. (2000) measurement, or through similar observations 
of other relaxed clusters. 

On the other hand, we note that dark halos of CDGs would steepen the central density 
profile and possibly complicate detection of a very small, fiat SIDM core. Ironically, the 
relaxed nature of rich clusters, which allows us to measure their masses by assuming hydro- 
static equilibrium, may turn out to obfuscate the details of the dark matter distribution. 
Even so, Chandra mass profiles appear to justify high-resolution N-body simulations with 
very low (cr < 0.1 cm^ g— 1) dark matter interaction cross sections. 

6. Summciry 

We have used a new spectral deprojection technique to derive a dark matter profile from 
a Chandra observation of EMSS 1358-1-6245. Our mass profile is consistent with optical weak 
lensing measurements (Hoekstra et al. 1998). It is nominally consistent with an NFW profile, 
although our derived concentration larger than expected from CDM simulations by a factor 
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of a few. Our best- fit NFW profile is characterized by = ISSlgg^ kpc and c = 8.4113. 
We also model the cluster as a non-singular isothermal sphere and place an upper limit on 
the core size Tc < 42 /i^q^ kpc (90% confidence). Comparing this value to the simulations of 
Yoshida et al. (2000), we rule out self-interacting dark matter with cross sections ^ 0.1 

2 —1 

cm g . 
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Fig. 1.— Broadband (0.3 < E/keV < 7.0) image of EMSS1358+6245. 




Fig. 2.— An HST image of EMSS 1358+6245 with superimposed soft X-ray (0.3 <E<2.0 
keV) contours. Each contour level is 3/5 of the previous value. 




Fig. 3.— Five USNO sources in the Chandra ACIS S2/S3 and HST fields. 
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Fig. 4. — Positional offsets of the three USNO sources which appear as point sources on the 
Chandra fields. 
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Fig. 5. — Broad-band X-ray contours overlaid on the adaptively smoothed Chandra X-ray 
image (top), the unsmoothed 2 x 2-binned image, and the HST field (bottom). 
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Fig. 6. — Wedge pairs used to estimate the cluster ellipticity. Wedges of equal radius are 
used first to determine the position angle, and then wedge B is scaled so that its luminosity 
is the same as wedge A. 
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Fig. 7. — Radial dependence of the flattening and position angle. It should be noted that 
these quantities are calculated by integrating to Ra, and in that sense represent a radially 
cumulative average. 



-28- 



r (h^o ^ kpc) 

50 100 200 500 800 50 100 200 500 800 



T 1 1 1 — I I I I 

15 - 



> 10 



5 



0.1 



I 



O 

I 

0) 

"XT 
"a 



0.01 



0.001 



T 1 1 — I — I — |— |— r 



(T)„=7.16±0.10 keV 



I ill 



li 



I I I I I 



J I I I I I I I 



1 1 — I — I — I — I I I 



T 1 1 1 — I — I — |— r 



0.1 0.2 0.5 1 



• • ••••• • 



_J I I I I I L_L_ 



2 0.1 0.2 
R (arc min) 



0.5 1 



Fig. 8. — Deprojected radial temperature and density profiles of EMSS 1358. Note that the 
inner two annuli are represented by two emission components at different temperatures. The 
profiles on the left are calculated allowing the hot gas temperature to vary; those on the 
right are frozen at 7 keV. 




Fig. 9. — Entropy profile of EMSS 1358. Tlie curve tliat spans tlie entire range in r represent 
tfie entropy of the gas at temperature Th; the short segment represents the gas at Tc. 
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Fig. 10. — Deprojected mass profile of EMSS 1358. The data points in the left panel 
represent M(r), the mass contained within spherical radius r, calculated using the density 
and unrestricted temperature profile shown in the left panels of Figure 8, with the single 
unphysical point represented as a dashed box. The solid dark blue curve and light blue 1-cr 
envelope represents the best-fit NFW profile. On the right the spherical NFW profile has 
been projected along the line of sight according to equation 15, in order to compare it to the 
weak lensing measurements of Hoekstra et al. (1998). 




Fig. 11. — The total integrated gravitating mass and gas mass of EMSS 1358. The non 
physical point (M(r) < 0) is represented by an open box in the total mass profile. 
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Fig. 12. — Power law fits to the mass profile of EMSS 1358. The top panel shows the mass 
profile for the = 7 keV model along with power law fits using points 1-k, k G [2,9]. 
The bottom panel shows the power law index of the fit, as a function of the radius of the 
outermost point used. Note that a is the power law slope on p(r), not M(r). 
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Fig. 13. — Softened isothermal sphere fits to a set of 5 simulated galaxy cluster observations 
which follow a P model profile. The fitted isothermal sphere core radius is plotted against 
the P model core radius. 
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Fig. 14. — Constraints on the SIDM cross section, adapted from Hennawi & Ostriker (2001). 
The hatch marks dehneate the region excluded by each constraint. The triangular region is 
derived from several astrophysical constraints of Hennawi & Ostriker (2001) (see text). The 
results of the simulations of Yoshida et al. (2000) and Dave et al. (2001) are indicated. The 
blue excluded region labeled 'Chandra limit' represents the constraint we derive from the 
absence of a significant core in the mass profile of EMSS 1358+6245. 
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Table 1. Source annuli statistics. 



annulus 






total 
photons 


source 
photons 


1 


11.5 


32.7 


2019 


2004 


2 


22.9 


65.2 


2004 


1958 


3 


34.7 


98.8 


2005 


1926 


4 


47.4 


135.0 


2007 


1895 


5 


61.6 


175.4 


2000 


1820 


6 


76.5 


217.8 


2000 


1761 


7 


93.0 


264.8 


2007 


1682 


8 


110.8 


315.5 


2011 


1589 


9 


153.8 


437.9 


4013 


2691 


10 


285.8 


813.8 


10040 


3299 
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Table 2. Comparison of hot and cold gas in the inner two 

shells. 



quantity component shell 1 shell 2 

T/keV cold 
hot 

p/10~^^ g cm~^ cold 

hot 

5M/W^ Mq cold 
hot 



1.83] 



-0.18 
-0.13 
7 07-1-5.40 
'■'5' -2.07 



-1-0.18 
0.13 



1.26 

7.0412-26 



19.4 



-0.5 



14.4 
2.57 



+1.5 
-7.2 
0.08 
0.06 



n. r-| O+0.013 
U.yzO_Q4g2 



n 264+° °^'^ 

o Qi +0.105 
•J-Ol _o.076 



V^/^total 



cold 
hot 



0.121 
0.879 



0.012 
0.988 



